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ABSTRACT 

The energetic balance of the Standard Solar Model (SSM) results from an equilib- 
rium between nuclear energy production, energy transfer, and photospheric emission. 
In this letter, we derive an order of magnitude of several % for the loss of energy in 
kinetic energy, magnetic energy, and X or UV radiation during the whole solar life- 
time from the observations of the present Sun. We also estimate the mass loss from 
the observations of young solar analogs which could reach up to 30 % of the current 
mass. We deduce new models of the present Sun, their associated neutrino fluxes, and 
their internal sound-speed profile. This approach sheds quantitative lights on the dis- 
agreement between the sound speed obtained by helioseismology and the sound speed 
derived from the SSM including the updated photospheric CNO abundances, based 
on recent observations. We conclude that about 20% of the present discrepancy could 
come from the incorrect description of the early phases of the Sun, its activity, its 
initial mass and mass-loss history. This study has obvious consequences on the solar 
system formation and the early evolution of the closest planets. 



Subject headings: Sun: helioseismology — Neutrinos — Sun: rotation — Sun: sur- 
face magnetism — Sun: UV radiation — Stars: Mass-loss 
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1. INTRODUCTION 



Impressive observational constraints unfolded during the past two decades to help check the 
physics of the solar interior. These effo rts led to a remarkable agreemen t between two complemen- 
tary p robes: neutrinos and seismology (|Turck-Chieze et al .1120041 12010 : iTurck-Chieze & Couvidat 
20101) . However, a reduction in the CNO solar photospheric composition h as now been confirmed 
by two independent groups (|Asplund et al.l 12009c ICaffau et al.1120081 120101) who directly extract it 
from photospheric and atmospheric lines, although a slight di f ferenc e exists between them. This 
long-expected revision (see section 2.1 of ITurck-Chieze et al.1 (| 1993b ) produced an inconsistency 
between Standard Solar Model (SSM ) predictions and the two aforementioned pro bes: neutrino 



fluxes and radial sound-s peed profile (ITurck-Chieze et al 



20081: iGuziketal. 



potheses were proposed (|Guzik et al.ll2005c iBasu & Antia 



2004 iBahcall et alJl2005h. Several hy 



2008c ITurck-Chieze. Nghiem & Mathis 



20 lOl) . Some of them have now been rejected, like a large error in the new pho- 



tospheric abundance estimate. Others are not currently easy to verify: 1) an incorrect microscopic 
diffusion: helium diffusion is, at least, measurable but the extraction of ra diation-zone CN O sig- 
natures in the equation of state casts some doubt on low subsurface values (|Lin et al.ll2007|) even if 
this information is more difficult to obtain; 2) role of a fossil magnetic field on the bound-bound 
contributions and of the atmospheric magnetic field on the photospheric composition; 3) underesti- 
ma ted opacity calculations: le ss than 10-15% per element if elements other than CNO are incorrect, 



see lTurck-Chieze et al.l (|201 1[) ... The latter is stimulating extensive opacity comparisons and spe- 
cific laboratory investigations (|Bailey et al.l 120071; ITurck-Chieze et al.ll2009l) . Unfortunately such 
experiments are not yet available for the whole solar radiative zone (RZ). 

In this letter, we examine another source of limitation of the SSM: the absence of solar ac- 
tivity effect during its evolution. Our paper on the transpo rt of momentum by rotatio n showed the 



limits of the SSM to predict the observed rotation p rofile (|Turck-Chieze et al. 



for including magnetic field in the stellar equations ( Duez. Mathis & Turck-Chiezel 2010 ). but 



20101) and the need 
we 



showed that their structural effects are small. Strong activity phenomena are now observed on 
young solar analogs, and consequently a transformation of nuclear energy into kinetic and mag- 
netic energy during the solar life must be studied. These phenomena produce a large variability 
in the X and UV spectra, inevitably accompanied by mass loss during the early phase of the solar 
evolution. The SSM does not take into account this dynamical aspect. Here we focus on the order 
of magnitude of such effects and on their past and present impacts. This letter is organized in three 
main sections. In section 2, we recall the SSM energetic balance and its limits. In section 3 we 
describe three updated solar models showing different types of loss of irradiance and mass loss 
connected to the past solar activity. The results are discussed in section 4. 
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2. THE ENERGETIC BALANCE OF THE STANDARD SOLAR MODEL 



At the end of the eighties, the pp n eutrino flux was easily estimated (ISpiro & Vignaudlll990h 
and the results of GALLEX and SAGE (lAnselmann et alll992[bavrin et al.lll992n could be inter- 
preted as a first evidence that neutrinos might oscillate during their travel from the solar center to 
the Earth. Indeed, a reasonable estimate of the pp neutrino flux reaching Earth is: 



O v pp = 2-^- * 



1 



^nucl 



And 2 



(1) 



where L nn r i is the energy produced by the pp chain, L is the s olar luminosity fixed at 3.846 x 
10 33 erg/s (|Guenther et al.lll992uTurck-Chieze & Couvidatll2010|) . and d is the Sun-Earth distance. 
This equation derives directly from the fact that: 



4p - > 4 He + 2e + + 2v e + 26.2 MeV 



(2) 



where 26.2 MeV is an approximate average value of the energy produced by the three chains 
ppl, ppll, and ppIII. From these equations, the pp neutrino flux reaching Earth is 6.49 10 10 cirT 2 /s. 
This estimate is correct at the 10% accuracy level, knowing that part of the solar neutrinos are 
also produced by the electronic capture of 1 Be and a small fraction (2 x 10" 4 ) comes from the 
disintegration of % B. The detailed estimate of the different neutrino fluxes was improved in the 
framework of the SSM which allows to solve the four e quations of stellar stru cture with extremely 



iramework ot me bbJvi wnicn allows to solve me lour e quations oi stellar stru cture witn extremel) 
refined physics. The Sun is evolved up to its present age (|Guenther et al. Il992l) of 4.60 + 0.04 Gyrs 



including 48 Myrs for the PMS. The solar luminosity increases by 30% during this evolution. The 
transport of energy by photons, throu gh interactions with the different species, is treated as properly 
as possible (|lglesias & Rogerslll996l) . Hydrostatic equilibrium is assumed during the evolution and 
no other source of energy is introduced. In the nineties, the predicted value of 6 x 10 10 cm~ 2 /s a t 
a 2% accuracy level was obtained by different groups, see table 13 of iTurck-Chieze et al.l (jl993|) . 
This predicted value remained nearly unchanged these last years. Therefore its only change results 
from the impact of the CNO cycle on the luminosity produced, which was slightly reduced by the 
recent CNO abundance update. 



Then, SNO definitely confirmed the existence of solar neutrino oscillations (|Ahmed et al 



20041) . measured a precise boron neutrino flux, and consequently provided an accurate determina- 
tion of the solar central temperature T c . The inclusion of neutrino oscillations and the CNO abun- 
dance revision led the SSM pp, 7 Be, and 8 1? neutrino predictions to be only in marg i nal ag reement 
with the neutrino detections (ICouvidat et al.l2003[ ITurck-Chieze et al.l2OOlaU2OQ4l.l2Ol 01) . On the 
contrary, the evolved Solar Seismic Model (SSeM), built to reproduce the observed sound speed, 
shows a remarkable agreement between its predictions and all the detection results (GALLEX- 
SAGE, CHLORINE, SNO, SUPER KAMI OKANDE and BOREXIN O), but it produces central 
solar conditions different from the SSM (see lTurck-Chieze et al.l(|2010|) and Table 1). 
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More and more, se ismic observations emph asize a dynamical view of the solar interior. The 



Hale cycle of 22 y ears (IHale & Nicholson 



1938) seems to result from motions of the whole con- 



vective zone (CZ) (|Brun. Miesch &Toomrel 120041: bikpati & Gilmanll2008l) . In MHD simulations 



of this zone, the radial luminosity at a fractional radius r, can be expressed as the sum of different 
components: 



L(r) 

47ZT 2 



= F e + F k + F r + F u +F v + F D 



(3) 



where F e ,F k ,F r ,F u ,F v ,F m are, re spectively, the enthalpy, kineti c, radiative, unresolved eddy, vis- 
cosity, and Pointing, fluxes (see (IBrun. Miesch &Toornrel 120041) for their definition). These terms 
are ignored in the equations governing the SSM. If most of them are believed to be small in the 
CZ, none has been estimated in the RZ, in particular for the young Sun that could rotate 10-100 
times faster than today, accompanied with related strong activity. 

Unfortunately, 3D simulations cannot yet produce an evolutionary model of the global Sun. 
Nevertheless, it is now well established that the present constant luminosity, name justified by 
the fact that in the SSM the solar luminosity var ies by only 10~ 8 during the last 100 years, is not 
confirmed by solar observations (|Frohlichl 120061) . The adopted present luminosity is constrained 
at the 10~ 4 level in our models. The observed variability is directly correlated with the solar 
Hale dynamo which produces a mean 10~ 3 luminosity variation, with an additional 3-4 xlO -3 
variation during the period of maximum activity of about 4-5 years due to the presence of preferred 
longitudinal spots appeari ng at the solar rotation periodicity. This variability can be accom panied 
by longer term variations (|Turck-Chieze & Lambertll2007t iTurck-Chieze & Lefebvrell201 1|) . 



3. A REVISED ENERGETIC BALANCE OF THE SOLAR INTERIOR 



The solar RZ encompasses ~ 98% of the total solar mass and its microscopic ingredients 
(nuclear reaction rates, opacity coefficients associated to the composition) imprint the sound-speed 
profile. The equilibrium between gravitational energy, nuclear energy production, and the energy 
escaping by photon interaction, is achieved in the radiative interior over long time scales. The 
sound-speed pro file down to 0.06 has been extracted from the Global Oscillatio ns at Low Fre 



quenc y (GOLF) (|Gabriel et al.lll995l) and the Michelson Doppler Imager (MP I) (IScherrer et al 



1995) instruments onboard SoHO (|Turck-Chieze et al 
cently been confirmed by BIS ON and MDI feasu et al 



200 lac ICouvidat et al.l 120031) . and has re 



2009j). In the RZ, a very slow meridional 



circulation could be in action (|Turck-Chieze et al.ll2010|) and a stable fossil magnetic field could 
still be present due to the very low diffusion coefficients considered. Both can affect the energetic 
balance, even though the introduction o f these processes in the stellar equat i ons show that their 



prese nt structural effects are really small (|Duez. Mathis & Turck-Chiezell2010l: ITurck-Chieze et al 



20101) . However a quantitative estimate of these dynamical processes on the energetic balance 
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requires to describe the early evolutionary stage, the evolution of the solar-core rotation, and the 
magnitude and topology of the fossil field. The present study attempts to give an order of magni- 
tude of these different effects. 

First, we estimate the maximum energy loss due to current radiative-zone dynamics. To this 
end, we use the central temperature difference between SSM and SSeM given in table 1. This cen- 
tral temperature is now determined with an accuracy of 5 x 10 3 thanks to the accuracy of neutrino 
fluxes. The pp luminosity varies like T 4 , where T c is the central solar temperature. Therefore, 
an increase in T c of 1.5% could be interpreted as an increase in central luminosity of 6% due to 
some energy loss by dynamical processes during the last million years, corresponding to the travel 
time of photons from solar center to the surface. This value is an upper limit because if the nuclear 
energy was partly transformed into another type of energy and evacuated, the Sun burnt more hy- 
drogen and its density was also slightly increased. Therefore we built two new solar models, for 
which we changed the calibrated luminosity at the present age by, respectively, 2.5% and 5%. 

The second estimate comes from the observation of young stars. In its in fancy, the Sun was 
certainly more active and emitted more strongly in X and XUV. In the review of iGudel I (|2007l) . the 
study of young solar analogs shows how the evolved X luminosity varies with time with respect 
to the total bolometric luminosity. This wavelength range represents between 5 x 10~ 7 an d 10~ 4 
of the total solar luminosity L© and varies today 10-20 times more. Moreover iRibasI (|20 1 0|) notes 
a variation of about 1000-2000 times the present emission in XUV for these young solar analogs. 
Therefore we applied the law described by these authors from their work on young stars: 



dh(f) 



ariGyrs) 



-1.23 



(4) 



where a is equal to 1.31 x 10 2 , and r is the age in Gyrs. This results in a luminosity loss of 
~ 3 x 10~ 2 at 0.5 Gyrs which is similar to the previous estimate when it is integrated over time 
because this loss decreases quickly with time. Thus the idea that the central solar region may be 
more evolved than what is today assumed in the classical stellar equations, is not such a stretch. 

Another consequence of solar activity is the associated mass loss. It is directly connected to 
the X luminosity. In order to see the consequences of this effect, we used the law suggested by 
Wood et al.l (|2005|) and we adopted the following decrease in mass for the first Gyrs: 



M w = 1.0 x 10 _1 l T(Gyrs) 



-2.23 



(5) 



as recommended by IRibasI (|2010|) . 



All the solar models of the p resent study include th e recent inputs on reaction rates (|Turck-Chieze et al 



2004J), the recent composition of lAsplund et al.l (|2009l) and the equation of state and opacity tables 
built with this composition. We also introduced the formalism suggested by lCanuto. Goldman & Mazitelli 
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Fig. 1. — Sound-speed profile in the solar core corresponding to the seismic sol ar model (solid black 
line) in perfect agreement with helioseismic data, to the standard solar model with lAsplund et all (120091) 
composition (solid red line), to a model with luminosity increased by 2.5% (dashed red line), and to a model 
with luminosity increased by 5% (dot-dashed red line). The vertical 3o~ error bars are very small and have 
been multiplied by 10 on this figure, while the horizontal error bars are rather large because few acoustic 
modes penetrate down to the core. 
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Fig. 2. — Squared sound-speed and density differences between the helioseismic inversions, obtained from 
the GOLF+MDI/SOHO acoustic modes, and various solar models: the standard solar model with Asplund 
(2009) composition (green double dashed line), the seismic solar model (black line) with the seismic error 
bars, and two models where a luminosity increase was introduced in the core (by 2.5 % in red dashed line 
and by 5 % in large dashed line). 
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Table 1: Boron neutrino flux prediction (expressed in 10 6 cm" 2 s _1 ) compared to the measured SNO 
value of 5.21 + 0.45, central temperature T c , initial helium content Y , central density p c , central 
composition in hydrogen X c and helium Y c , surface helium Y s , and a parameter, for the seismic 
solar model, and the standard and non-standard solar models. All the SSM are computed with the 
Asplund (2009) composition. 



Model 



Boron neutrino flux 



Pc 



SSeM 



5.31+0.6 
X c = 0.339 



15.74 
Y c = 0. 645 



0.277 
a MLT = 2.04 



153.02 
Y,= 0.251 



SSM 



4.50 15.54 0.2645 150.6 

X c = 0.357 Y c = 0.627 a CGM = 0.786 Y s = 0.235 



SSM L+0.025 



5.32 15.67 0.2670 153.8 

X c = 0.347 Y c = 0.637 a CGM = 0.810 Y s = 0.237 



SSM L+0.050 



6.33 15.79 0.269 157.0 

X c = 0.338 Y c = 0.647 a CGM = 0.836 y s = 0.239 



SSM Tac 



4.2 15.49 0.2671 150.3 

X c = 0.362 y c = 0.623 « C gm= 0.712 ^=0.240 



SSM M ;mf =1.33 



4.74 15.58 0.2603 154.1 

X c = 0.347 y c =0.637 a CGM = 0.829 y s = 0.234 
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(I1996Q that improves the description of convection (see references in IPiau et al.l (| 20111) ). All the 
results are summarized in Table 1 and Figures 1, 2, and 3. 

Table 1 lists the models computed. The evolved SSeM, built to reproduce the observed 
sound speed, agrees with all the observables: neutrino fluxes and photospheric helium abundance. 
The new SSM shows a clear disagreement for the photospheric helium content. The discrep- 
ancy on helium is reduced by 20% with the introduction of mixing in the tachocline following 
Brun. Turck-Chieze & Zahnl (|1999|) . The SSM models with a central luminosity increased by, 
respectively, 2.5 and 5 % through a change in calibration show a net increase in the central temper- 
ature, and consequently an improvement of the boron neutrino predictions even though the second 
case overestimates the neutrino flux. The mass-loss impact is less visible on the temperature than 
on the composition. 



Figure 1 shows the so und speed in the c ore obtained from GOLF (|Turck-Chieze et al.ll2001al) 
and confirmed by BiSON (|Basu et al.ll2009|) . While the vertical error bars are extremely small, 
the horizontal ones remain large in the region below 0.3 R Q which contains more than half of the 
solar mass. This uncertainty could be red uced by measuring mixed-mode frequencies of different 
degrees or high-frequency gravity modes dMathur et al.ll2007l) . By construction the seismic model 
(full black line) is in perfect agreement with the inverted sound speed from helioseismic data. We 
added three solar models: the SSM with turbulence in the tachocline and the two models assuming 
a loss of luminosity of 2.5 and 5% during the solar lifetime. These models (as shown in table 1) are 
in better agreement than the SSM in the actual core where the boron neutrinos are emitted. This 
figure shows the interest to better describe the solar core. More constraints in the inner part are 
needed not only to follow the dynam ical processes but also to see if some dark-matter signature 
can be detected (ILopes & Silldl2010h which will differ from the signatures that we look for in the 
present study. 

Figure 2 shows the squared sound-speed and density differences between helioseismic values 
and various solar models for the entire solar interior. The differences for the SSM (dotted line) 
using the new photospheric composition and the full spectrum of turbulence model of Canuto et al. 
(1996), remain important. In order to separate the effect of luminosity loss from the effect of mass 
loss, we show the models where the luminosity is increased by +2.5% in small dashed line and by 
+5% in large dashed line. The effect is not large and has a sign opposite to the one in the core (see 
Figure 1). Using Equation 4 for the luminosity loss at the surface during the early stage, leads to 
a smaller effect because the star readjusts very quickly at each loss. The comparison of the two 
approaches shows that during the solar life an energy loss of a few % is possible. However these 
simulations are too crude as the energy produced is not explicitly transformed into another energy 
type. Therefore Equation 4 alone cannot reproduce the actual RZ. 



In addition to SSM (in large dashed blue) and SSeM (in black), Figure 3 recalls the influence 
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Fig. 3. — Same as Figure 2. The updated standard solar model is in green double dashed line, the SSM 
model with mixing in the tachocline is in red dashed line, the seismic solar model is in black line with the 
seismic error bars, and a model with mass loss in the early evolutionary phases following observations from 
young solar analogs is in blue long dashed line. 
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of the introduction of mixing in the tachocline ( Brun. Turck-Chieze & Zahnlll999 : Piau & Turck-Chieze 



2002|) . The horizontal turbulence in the tachocline slightly shifts the peak discrepancy downwards 



and suppresses the peak in density at the transition between radiation and convection. It reduces 
the discrepancy in the photospheric helium by 20% (see Table 1), but slightly decreases the cen- 
tral temperature. These points have already been mentioned in our previous works. Figure 3 also 
shows the influence of mass loss given by Equation (5) on the present squared sound-speed and 
density differences. This mass loss is applied starting at 50 Myrs, corresponding to the beginning 
of the main sequence, and mimics the result of magnetic wind. The initial solar mass is 1.33 M Q 
(where M Q is the present solar mass). An even higher initial mass could be considered if the mass 
loss was applied at the separation of the solar disk near 3 Myrs. In the present study, the luminosity 
at the beginning of the main sequence is 1 .5 L Q instead of 0.7 L Q , and the young Sun reaches 1 M Q 
at 0.77 Gyr. Its luminosity is smaller by 10% than the present luminosity at 1.5 Gyrs instead of 
20% with the SSM. In that case, the discrepancy on the squared sound-speed and density profiles 
is reduced and the predicted boron neutrino flux is slightly increased (table 1). 



4. SUMMARY, DISCUSSION AND CONSEQUENCES 

This paper confirms the difficulty for the SSM to appropriately reproduce all the present solar 
observations: neutrino fluxes, sound-speed and density profiles, photospheric helium content, and 
location of the base of the convective zone. This paper proposes two directions of investigation 
related to the activity of the young Sun and to the comparison of central temperature between SSM 
and SSeM. We derive the order of magnitude of the corresponding effects and the quantitative 
implications on the aforementioned observables: (1) the Sun could have transformed about 2.5-4 
% of the energy produced during the early evolutionary phases into another form of energy due to 
strong activity and the related interior dynamics, with a redistribution (or evacuation) of nuclear 
energy through kinetic and magnetic energies in the RZ; and (2) its initial mass could have been 
greater by 20-30% than the present solar mass if the relation observed on young solar analogs is 
used. Then, some residual effects on the present observables cannot be ruled out if we compare the 
difference of squared sound speed, the central temperature between SSM and SSeM and the boron 
neutrino fluxes. 



These ideas are not new (ITurck-Chieze et al.lll988t, ISackman & Boothroydll2003l : iGuzik et al. 



20101) but, for the first time, they are deduced from observations of young solar analogs and from 
the comparison of two solar probes, and are not only parametrized. Of course the mass loss studied 
here is rather large in comparison with previous works (limited by the previous agreement on the 
sound speed and lithium depletion) but is not ruled out by observation of young solar analogs. 
The knowledge of the present rotation rate in the solar core could help put more constraints on 
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the phenomena studied here. This study will be pursued by a de dicated work on the histo ry of 
the internal magnetism after introduction of the rotation history (|Turck-Chieze et al.ll2010|) . An 
active young Sun probably possesses a convective core and an early internal dynamo in addition 
to some mass loss: this would greatly modify the initial condit ions of planet formation and could 
partly solve the "solar paradox" mentioned a long time ago by ISagan & Mullerl (|1972|) and never 
really t otally solved: how such a low initial solar lumi nosity could be compatible with Mars early 
history (|Forget & Pierrehumbertlll997[ IHaberle Ill998|) ? This young Sun an d solar analogs would 



also a ffect the poorly understood history of 7 Li during their early evolution (IPiau & Turck-Chieze 
2002k 



Moreover, other ideas must be explored like the magnetic field effec t on solar spectral li nes 
which may impact the composition determination in the solar atmosphere (|Fabbian et al.ll2010i) . or 
the effect of a potential fossil field on opacities in the RZ, or the role of gravity waves. 



We thank the anonymous referee for his (her) suggestions and improvements of the first ver- 
sion. 
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